Voltage profile and four terminal resistance of an interacting quantum wire. 
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We investigate the behavior of the four-terminal resistance Ript in a quantum wire described by 
a Luttinger liquid in two relevant situations: (i) in the presence of a single impurity within the wire 
and (ii) under the effect of asymmetries introduced by disordered voltage probes. In the first case, 
interactions leave a signature in a power law behavior of Ri p t as a function of the voltage V and 
the temperature T. In the second case interactions tend to mask the effect of the asymmetries. In 
both scenarios the occurrence of negative values of R^ v t is explained in simple terms. 

PACS numbers: 72.10.-Bg,73.23.-b,73.63.Nm, 73.63.Fg 
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The proposal of a fundamental relation between the 
two terminal conductance G and the universal quantum 
Go = e 2 /h is one of the milestones of electronic quantum 
transport in mesoscopic systems^. For non-interacting 
electrons, such a relation explicitly reads G = uGq, being 
n the number of electronic channels^. The same relation 
has been later theoretically and experimentally probed to 
be also valid in the case of wires of interacting electrons of 
finite length ideally attached to non-interacting leads^. 

Experiments in single- wall nanotubes (SWNTs), 
ropes of SWNTs and also in multiple-wall nanotubes 
(MWNTs)^ have, instead, identified that the tunneling 
conductance to metallic contacts, G t follows a power law 
behavior with the voltage V G t oc V a at low temperature 
T, and G t oc T a for low V. The exponent a being a func- 
tion of the forward electron-electron (e-e) interaction g. 
These features can be understood within the framework 
of a Luttinger liquid (LL) theory by means of theoretical 
treatments^ going beyond linear response in V and T. 

Recently, a combined structure of MWNTs and 
SWNTs has been used to analyze the behavior of the 
four-point resistance Ri p t of a SWNT 8 . The total re- 
sistance of a mesoscopic system in a two terminal setup 
contains the component X/Gq due to the coupling to the 
reservoirs. Instead, R^ v t is expected to characterize the 
genuine resistance of the sample. For non-interacting 
electrons at low V and zero temperature, Biittiker has 
elaborated the concept of the multi-terminal resistance 
within scattering-matrix theory (SMT) 7 , emphasizing 
the role of the symmetries. Although a naive expectation 
would be i?4 pt > 0, this theory predicts also the possibil- 
ity of Ri vt < as a consequence of quantum interference 
effects. This remarkable feature has been experimentally 
observed^. 

While the consequences of elastic scattering due to im- 
purities can be analyzed in terms of non-interacting elec- 
trons, the role of the e-e interaction in the behavior of 
Ript remains an open question. The proper evaluation of 
this quantity implies dealing with a multi-terminal setup 
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FIG. 1: Sketch of the setup: A voltage V is imposed on a 
Luttinger liquid with a backscattering impurity of strength 
As, through the chemical potentials for the left and right 
movers: jj, Tt i = /i± V/2. Two voltage probes are connected at 
the positions xi,Xa. The corresponding chemical potentials 
fii t 2 are fixed by the condition of zero current through the 
contacts. 



as the one of Fig. [TJ which is difficult to implement within 
theoretical approaches like those of Refs. m MTolfTl] . Previ- 
ous multi-terminal treatments in LL rely in the effective 
reduction to a non-interacting model by recourse to a 
Hartree-Fock decoupling of the interaction termi 2 . or fo- 
cus in linear response in V—. In this work we use Keldysh 
non-equilibrium Green's functions, which is a convenient 
framework to tackle multi-terminal geometries by exactly 
treating the e-e interaction while going beyond linear re- 
sponse. We analyze two relevant ingredients giving rise 
to a non-trivial behavior of the local potential and 
Ript- (i) the presence of an impurity in the wire and (ii) 
a clean wire probed by disordered leads with asymmetric 
densities of states. 

We consider the setup sketched in Fig. [TJ which 
consists in a quantum wire described by a Tomonaga- 
Luttingcr model under the influence of a voltage V. Two 
reservoirs are in contact with the wire at the points x\ , 
through very weak ("non invasive") tunneling constants, 
their chemical potentials /xi, /i2 satisfy the condition of a 
vanishing current through the ensuing contacts. In this 
way, a continuous current I flows through the wire, and 
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In the case of an ideal and clean wire with identical non 



2 



invasive probes, a simple analysis of the symmetries of 
the setup leads to R± pt = 0. The full system is described 
by the action S = S W i re + Si m p + S res + S con t-, where S res 
describes the two reservoirs that constitute the voltage 
probes and S W i re is expressed in terms of right (r) and 
left (I) movers (in our system of units % = Vf = e = 1): 

Swire = J dx dt {lpl{i(d t + d x ) - Urtyr + 

+$[i(fit ~ d x ) - m]iln - g[4^ r + V^/] 2 },(2) 

where g is the Luttinger e-e interaction in the forward 
channel while p, r = p + V/2 and pi = p — V/2. The 
effect of the impurity is contained in the back-scattering 
interaction: 



S inp = A B / dxdt5{x-b)[e- 2 ^ x ^% + H.c].{?>) 



The term S cont represents the tunneling between the 
reservoirs and the wire: 



Scont ^ 

j=l,2,a,P=l,r 



dxdtwjS(x — Xj) x 



H.c] 
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where the fields \ a j> with a, f3 = I, r, j = 1, 2 denote de- 
grees of freedom of the reservoirs. The upper and lower 
sign corresponds to a = r, I, respectively. For simplicity, 
the dependence of the fields on x and t has been omitted 
in the above equations. At this point we carry out the 
gauge transformation ipj r (x) — > e i±kpx ip a (x) where fcp 
is the Fermi vector of the electrons in the wire (a simi- 
lar transformation involving k' F is implemented with the 

fields Xaj)- 

The tunneling currents through the contacts to the 
probes read 
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In what follows, we evaluate the currents Ij up to the first 
order of perturbation theory in the tunneling amplitudes. 
This procedure is appropriate in the limit of weak Wj, 
which is a reasonable assumption for measurements of 
i?4 P t with 'non-invasive' probe leads&. Within this lowest 
order of perturbation theory: 
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-G<^{x hXj] u>)Gf{u>)], (6) 



where G^'^(u>) are the Fourier transforms with respect 
to t — t' of the lesser and bigger Green's functions 
G< p {x,x';t-t') = i(ipl(x',t')ip a (x,t)), G> (x,x';t- 
t') = —i(il; a (x,t)ipp(x',t')), corresponding to the 
wire uncoupled from the probes, while G - ' (a;) = 



Aj >,< (w)pj(a;) are the Green's functions of the uncou- 
pled probe reservoirs, with Xf(uj) = if(u> — fXj), \f(w) = 
— f (uj — Pj)], and pj(u>) the ensuing density of states, 
which we assume to be identical for the two kinds of 
movers within these systems. The chemical potentials of 
the probes, p,j must be fixed to satisfy Ij = 0. Notice 
that, within this order of perturbation theory, the effect 
of the two probes is completely uncorrelated from one 
another, since interference terms between the probes and 
resistive effects involve second order processes in wj^. 

We now turn to analyze the first ingredient of interest, 
namely, the effect of an impurity in the wire. In terms 
of our model, this corresponds to consider a finite As. 
We also consider a simple model for the probes, with a 
constant density of states pj (us) = po . While the Green's 
functions of the uncoupled homogeneous interacting wire 
are knowni^, the evaluation of the corresponding func- 
tions in the presence of a backscattering center is a non- 
trivial task. Below, we indicate the lines we have followed 
in order to evaluate them up to the first order of pertur- 
bation theory in As. The expressions for the lesser and 
bigger Green's functions cast: 



G>'<(x,x',lj) 



S a p^a' < ( U )p0,a{ X ~ x',Uj) + 

<5a/jA. B {A>' < (w)po,a(z - b,uj) 
x[G^-b,u)]* + X>'<{u) 
xGX(x-b,uj)p o , (b-x',cj)},(7) 



with I = r,r = I, A<(w) = if(u> - p a ), A>(w) = -i[l - 
f{uj — p a )]. The spectral density 

p , a (x,uj + p a ) = C,p exp [=R(— - k F )x] 

v 

xM 2 X 7 ,27+l,±2^-), (8) 
v 

corresponds to the clean LL uncoupled from the probes, 
where v = \J\ + 2g/ir is the renormalized Fermi velocity 
and 4>(a, b; c) is Kummer's hypergeometric function. The 
exponent 7 = (K + K^ 1 - 2)/4, (K = 1/v) is deter- 
mined by g. The retarded Green's functions are defined 
from the Kramers-Kronig relation, being po t0l (x,u)) = 



-2Im[G« (3 



In order to perform numerical com- 



putations we introduce an energy cutoff A by replacing 
po,a(x, u + p a ) — > 0(|w| — A)po ia (x, u> + p a ). Therefore 
the constant is a function of A which is determined 

by the sum rule J^°° du!p a fi(0,uj) = 2tt. 

Substituting the above expressions in © gives the fol- 
lowing result for the currents through the contacts: 

/DO 
— [f(uj- p-j) -f{w-fi a )] 

x Pj (u)plff{ Xj ,Lj), (9) 

being 

p e J f ( Xj ,oj) = p , a (0,u) +2\ B Re{po !a (xj - b,u) 

xK,n(xi-b,")]*h (10) 
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FIG. 2: Behavior of the local potential fjbj for a wire with 
As = 0.1 as a function of the bias V for T = and Xj = —0.2 
(upper left panel) and as a function of T for V = 0.2 and 
Xj = —0.2 (right panel). Circles, squares and triangles corre- 
spond to K = 1, 0.5, 0.4, respectively. Solid and open symbols 
correspond to V = 0.05,0.2, respectively (right panel). The 
fits with power laws are shown in dotted lines. Lower panel: 
fij as a function of x — b for V = 0.2, . . . , 1 and K = 0.5. 
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FIG. 3: Four point resistance K^ v tjRi v t of a LL with interac- 
tion parameter K for a bias p T — = V . The probe j ' = 1 has 
a Lorentzian density of states (a resonance) with width A = 
0.4V centered at Eq. The probe j = 2 has a constant density 
of states (which implies fi% — fj, = (/i; +/i r )/2). The left panel 
corresponds to temperature T — and shows the changes in 
Ript as the center of the resonance is moved around p. Cir- 
cles, squares and triangles correspond to K = 0.4, 0.5, 0.6, 
while dashed lines, to the non-interacting case (K = 1). The 
right panel shows the evolution with the temperature for the 
value of (Eo — p) /V = —0.3 indicated with dotted lines in the 
left panel. The thin dashed line corresponds to K = 1 and 
{E - p)/V = 0.3. 



the effective density of states of the a movers (or the 
a-injectivity-ii) at the position Xj of the wire, which con- 
tains the contribution of the local density of states of the 
homogeneous wire, po,a(0, to) plus a correction due to the 
backscattering by the impurity. It is important to note 
that Ij(xj) oc Xb, while the current through the wire is 
I oc A| (see also Refs^SJi). For T = and low V, the 
condition /, = leads to 



V pi ff (xj,fi) - p%H{x$,\i) 



^3 >= M + 



efft 



2 Pi 



Pr ff (xj,p,) 



(11) 



where we have made use of the fact that pj (u>) is constant 
and that p e Jf (xj , to) ~ p^* (xj , p) within a small window 
Jw - p\ < V/2. Eq. CP reduces to Eq.(37) of Ref. 
[IT] for non-interacting electrons. It explicitly shows that 
the local potential monitors the difference between the 
left and right injcctivities at the point Xj, relative to the 
total density of states at the given point. It is natural to 
expect that such an observable should provide valuable 
information on the Friedel oscillations introduced by the 



impurity 
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, as well as on the strength of the e- 



e interactions. In fact, a low-energy expansion of the 
spectral densities casts: 

Pj - m + Ci As sm(2k F x)V 2 "« + \ T = 

Pj - p + C 2 X B Vsm{2k F x)T 2 '< 1 V«r, (12) 

with Ci, Ci functions of 7. 

The results of the full numerical calculation of pj from 
the condition of Ij(xj) = for arbitrary voltage differ- 
ences and temperature, with the exact effective density 



p%ff(xj,u>) (also evaluated numerically from eq. (fT0|) ). are 
shown in Fig. El At T = 0, function of the dis- 

tance to the impurity, shown in the lower panel, oscillates 
with the period 2k F of the Friedel oscillations, with a 
voltage-dependent amplitude that follows the power law 
(fT2"]l (see upper left panel). For large distances to the im- 
purity and high V , beyond the scope of the approxima- 
tions leading to eqs. (TH?)) . the pattern shows additional 
structure, and the amplitude of the oscillations decreases 
with the distance to the impurity. The evolution of pj as 
the temperature increases, corresponding to two values 
of the voltage, is illustrated in the upper right panel of 
Fig. [21 where the dependence (JTSJ) is also verified within 
the low T and low V regime, with V -C T. From these 
features we can infer the behavior of i?4 P t /i?2pt along the 
sample by simply substituting (lll| in |T]). In particular, 
we conclude that, as a function of x, Ri p t/R2pt should 
follow the pattern of Friedel oscillations, being positive or 
negative, depending on the points at which the probes are 
connected. As a function of V it should be a power law 
with exponent 2j. As a function of T, it should present 
rapid changes within the range T < V and a crossover to 
a power law with exponent 2j at higher temperatures. 

We now consider the second situation of interest: the 
wire without impurities (Xb =0) but disordered probe 
leads. Thus, the expressions for the Green's function of 
the wire reduce to the ones for the homogeneous LL while 
the expression for the tunneling current of eq. © reduces 
to eq. (HI) with p e Jf(xj,w) = p a ^(0,Lu). While perfect 
metallic systems are expected to have approximately flat 
densities of states, impurities introduce effective barriers, 
generating peaks in the densities of states pj(u>)- 
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Let us first analyze T = and probes with asymmetric 
density of states, such that Pj(u) ~ pj(p),p — V/2 < 

u> < p and Pj(i^j) ~ Pj(fi), p < w < /i + V/2. Then, the 
condition of a vanishing current © leads to: 



Pj = P 



v p+(pY -p-(pY 



(13) 



being v = 1/(27 + !)• F° r probes with symmetric den- 
sities of states, we get /«j = p and i?4 P t = 0. Instead, 
for an asymmetric density of states, the potential drop 
between the highest potential r and the probe j is lower 
(higher) than the one between j and I for pj (p) > pf (p) 

(pj(p) < Pj(p)), respectively, which reflects the fact 
that the larger the spectral weight of the probe, the larger 
the ability of that element to introduce resistive effects. 
For finite temperature and very low voltage such that 
T 3> V, it can be verified that pj — p,j — 1,2 and 
Ript = 0. 

Therefore, asymmetric densities of states of at least one 
of the probes together with the condition pi(u>) =fc Piiu) 
would lead to a non vanishing i?4 P t when T < V . An 
example is analyzed in Fig. [3] We consider a Breit- 
Wigner model for one of the probes, assuming a single 
resonance within the window of width V centered around 
p: pi(ui) — A\/[{uj — E ) 2 + A], and a constant density 
of states p2 — A2 for the other probe, where A\ , A2 are 
normalization constants and p — V/2 < Eq < p + V/2. 
Under these conditions P2 = p, while p\ is determined 
to satisfy 1\ = 0. Results for the corresponding relative 
resistance Ri p t/R2 P t — {pi — P2)/V are shown in Fig. 
[3J The left panel corresponds to temperature T = 0. 
When the center of the resonance Eq coincides with the 
mean chemical potential p, the spectral weight spreads 
out symmetrically around this point. Thus, Pi(p) = 
Pi(p) and pi = p, then Ri p t — 0. As the center of 
the resonance moves to lower energies, so does p\ and 
R^pt becomes negative. Conversely, for Eq > p, it is 
obtained Ri p t > 0. Remarkably, interactions tend to 



mask the structure observed in the non-interacting case 
(with K — 1). The behavior of Ri p t as a function of the 
temperature is shown in the right panel for the case (Eq — 
p)/V = -0.3. Notice that the cases with (E -p)/V > 
can be obtained from the ones with (Eq — p)/V < by 



—Ript in the figure. In all 



simply transforming Ri p t 
the cases, there is a range of temperature T < V, where 
Ript experiments significant changes. 

To conclude, let us comment on the theoretical and 
experimental impact of our results. When an impurity is 
in the wire, it induces Friedel oscillations that manifest 
themselves in the local voltage and Ript- The interac- 
tions leave a clear signature in the power law behavior 
Ri P t/R2pt oc V 21 and Ri pt /R2 P t °c T 27 . Interestingly, 
the exponent is different from the one predicted in Ref. 
[ill for the two terminal conductance of a LL with an 
impurity. This result has a significant conceptual weight 
since it constitutes a concrete example of the fact that dif- 
ferent fundamental processes contribute to each of these 
quantities (Ri p t/R2pt k As while G oc \ 2 B ). Therefore, a 
genuine multiterminal setup is essential to evaluate i?4 pt . 
For impurities in the probes and an asymmetric configu- 
ration, i?4 P t is determined by the way in which the den- 
sity of states of the probe is distributed within an energy 
window of width V centered in p, while the e-e inter- 
actions play a milder role. In both cases, the behavior 
of Ript as a function of temperature at a sizable V , is 
highly non universal and exhibits significant changes in 
the range T < V . These results should help to provide 
a theoretical framework to further analyze experimental 
data in SWNT, like those of Ref. |8| as well as to guide 
additional experiments along that line in the future. 
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